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ABSTRACT 

' We study alignment of grains by radiative torques. We found steep rise of radiative torque efficiency 

as grain size increases. This allows larger grains that are known to exist within molecular clouds 
. to be aligned by the attenuated and reddened interstellar radiation field. In particular, we found 

that, even deep inside giant molecular clouds, e.g. at optical depths corresponding to Ay < 10, 
large grains can still be aligned by radiative torque. This means that, contrary to earlier claims, 
far-infrared/submillimctcr polarimetry provides a reliable tool to study magnetic fields of pre-stellar 
cores. Our results show that the grain size distribution is important for determining the relation 
between the degree of polarization and intensity, 
^vq ' Subject headings: ISM: dust, extinction — ISM: clouds — polarization radiative transfer 

1. INTRODUCTION 

' It is widely believed that magnetic fields play a crucial role for the dynamics of molecular clouds and for the star 
ly-j , formation processes (see review by Crutcher 2004 and references therein). 1 One of the most informative techniques of 

■ studying magnetic fields in molecular clouds is based on the use of starlight polarization and polarized emission arising 
from aligned dust. 

Alignment of interstellar dust was not expected by theorists. Very soon after the discovery of the interstellar origin 

■ of starlight polarization by Hall (1949) and Hiltner (1949), it became clear that interstellar grains get aligned with 
respect to magnetic field. It did not take long time to realize that grains tend to be aligned with their long axes 
perpendicular to magnetic field. However, progress in theoretical understanding of the alignment has been surprisingly 

I ■ slow in spite of the fact that great minds like L. Spitzer and E. Purcell worked on grain alignment (see Spitzer & Tukey 
1951, Purcell 1969, Spitzer & Purcell 1971, Purcell 1975, 1979, Spitzer & McGlyn 1979). The problem happened to be 

+3 , very tough and a lot of relevant physics had to be uncovered. An extended discussion of different proposed mechanisms 

^ ■ with the relevant references can be found in a recent review by Lazarian (2003). 

Originally it was widely believed that interstellar grains can be well aligned by a paramagnetic mechanism (Davis 
. & Greenstein 1951). This mechanism, based on the direct interaction of rotating grains with the interstellar magnetic 
field, required magnetic fields that are substantially stronger than those uncovered by other techniques. 2 Later, a 

r> | pioneering work by Purcell (1979, henceforth P79) showed a way how to make grain alignment more efficient. Purcell 
noticed that grains rotating at high rates are less susceptible to the randomization induced by gaseous collisions, 
while paramagnetic alignment would proceed at the same rate. He introduced several processes that are bound to 
make grains very fast " suprathermal" rotators. They are a) variations of the accommodation coefficient for atoms and 
molecules bouncing from the grain surface, b) variations of the coefficient of electron ejection, c) variations of the sites 
of H2 formation over grain surface. As H2 formation over grain surfaces is a common interstellar process and every 
formation event could deposit an appreciable angular momentum with the grain, Purcell identified c) as the major 
cause of grain fast rotation. He claimed that the catalytic sites ejecting H2 molecules ( frequently called "Purcell's 
rockets") should spin up grains very efficiently for most of the diffuse ISM. One can easily see that within the Purcell's 
model, even a small fraction of atomic hydrogen present in molecular clouds would also make them suprathermal (i.e. 
Ekin 3> kTg ra i n ). For decades this became a standard explanation for grain alignment puzzle, although it could not 
explain several observational facts, e.g. why observations indicate that small grains are less aligned than the large 
ones. 

New physics of grain internal motion uncovered fairly recently explains why small grains are poorly aligned by 
Purcell's mechanism. The inefficiency stems from grain internal wobbling. Indeed, for sufficiently small grains it is 
impossible to assume that they rotate perfectly about their axis of maximal inertia. It is interesting to recall that 
the issue of grain wobbling was a part of the alignment process discussed e.g. by Spitzer (see Jones & Spitzer 1968). 

1 Existing claims of the contrary (see Padoan & Norlund 1998) make quantitative studies of magnetic fields more essential. 

2 As discussed, for instance in Lazarian (2003), the very small grains are likely to be aligned by this mechanism and this can explain the 
peculiarities of the UV part of the spectrum of the polarized radiation observed (see Kim & Martin 1995). The alignment of small grains 
by paramagnetic relaxation is possible as efficiency of the Davis-Greenstein mechanism increases as the grain size decreases. The degree 
of alignment of small grains provides a direct constrain on the intensity of magnetic field, which is a subject that calls for more of UV 
polarimetry work. 
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However, when P79 identified Barnett relaxation as an fast process of internal relaxation that aligns grain rotation 
with the axis of the maximal moment of inertia 3 , an idea that grains always rotate about the axis of maximal inertia 
got universally accepted. The flaw with this reasoning was found in Lazarian (1994), where it was shown that thermal 
fluctuations within grain material induce grain wobbling, the amplitude of which depends on the ratio of the grain 
rotational energy E/-in and kT gra i n . 4 The quantitative theory of the effect presented in Lazarian & Roberge (1997) 
allowed to revise the Spitzer & McGlynn (1979) theory of crossovers (Lazarian & Draine 1997) as well as the theory 
of paramagnetic alignment of thermally rotating grains (Lazarian 1997a, Roberge & Lazarian 1999). 

However, a more interesting development came about later, when Lazarian & Draine (1999a, henceforth LD99a) 
realised that grains not only wobble but occasionally flip. LD99a found that small grains flip more frequently that 
large ones. As the result regular torques, e.g. torques due to ejection of H 2 molecules, get averaged out over flipping 
grains and they get "thermally trapped", i.e. rotate at thermal velocities in spite of the presence of Purcell's torques. 
Taking into account that the paramagnetic alignment of thermally rotating grains is rather inefficient (see Roberge & 
Lazarian 1999) it is possible became possible to explain why small grain may be poorly aligned. 

A new twist to the theory of grain alignment came about when Lazarian & Draine (1999b, henceforth LD99b) found 
that species with nuclear moments within the grain, e.g. 1 H, 13 C, 27 A1, 31 P, 29 Si, Mn... , bring about a new type 
of internal relaxation which was termed "nuclear relaxation". This type of relaxation for grains larger than 10~ 5 cm 
happened to be ~ 10 6 times more efficient that the Barnett relaxation introduced in Purcell (1979). As the result, 

LD99b claimed that for diffuse interstellar gas nuclear relaxation thermally traps grains of the sizes from 10 cm 

making the Purcell mechanism inefficient. 

A group of alternative mechanisms of alignment that rely on the relative gas-grain motion have their particular 
niches. The first mechanical alignment mechanism was pioneered by Gold (1951). Later work included driving grains 
by ambipolar diffusion (Roberge & Hanany 1990, Roberge, Hanany & Messinger 1995) and Alfvcn waves (Lazarian 
1994, Lazarian 1997b, Lazarian & Yan 2002, Yan & Lazarian 2003). Although new efficient processes of mechanical 
alignment have been proposed (Lazarian 1995, Lazarian & Efroimsky 1996, Lazarian, Efroimsky & Ozik 1996), this 
did not make mechanical alignment universally applicable. 

All this provided the background that made radiative torques mechanism most promising for explaining grain 
alignment over vast expanses of the interstellar space. Introduced first by Dolginov (1972) and Dolginov & Mytrophanov 
(1976) the radiative torques were mostly forgotten till the pioneering work by Draine & Wcingartncr (1996; hereafter 
DW96), where their efficiency was demonstrated using numerical simulations. The radiative torques make use of 
interaction of radiation with a grain to spin it up. Indeed, in general, one would expect that the cross sections of the 
interaction of an irregular grain with left and right circular polarized photons are different. As the non-polarized light 
can be presented as a superposition of the equal fluxes of photons with opposite circular polarization, the interaction 
of such a light with the irregular grain would result in grain spin up. Unlike Purcell's torques, that are fixed in grain 
frame, the radiative torques are expected to be less affected by grain flipping. 

The predictions of radiative torque mechanism are roughly consistent with the molecular cloud extinction and 
emission polarimetry (Lazarian, Goodman & Myers 1997) and the polarization spectrum measured (see Hildebrand et 
al. 2000). They have been demonstrated to be efficient in a laboratory setup (Abbas et al. 2004). Evidence in favor 
of radiative torque alignment was found for the Whittet et al. (2001) data obtained at the interface of the dense and 
diffuse gas at the Taurus cloud (see Lazarian 2003). 

In view of this success the radiative torque mechanism is the primary mechanism that we are going to study in 
relation to grain alignment deep inside molecular clouds. While a possible failuc of radiative torques there does not 
exclude that grains arc aligned deep within molecular clouds, their success would definitely make polarimetric studies of 
molecular clouds much more trustworthy and informative. Whether grains are aligned there is necessary to understand 
to know whether aligned grains trace only surface magnetic fields or magnetic fields deeply embedded into molecular 
clouds. The topology of magnetic field inside molecular clouds is essential for understanding for star formation. 

It has been shown that optical and near infrared polarimetry provide magnetic fields only to A v of 2 or less (Goodman 
et al. 1995, Acre et al. 1998). Is it the same for far infrared polarimetry? This is the question that we address in this 
paper. Earlier answers (see Lazarian et al. 1997) appeal to stars embedded in the cloud. Indeed, such stars can induce 
alignment through their radiation. Here we consider an extreme case of a cloud without any embedded stars. This 
situation is also motivated observationally, as some recent observations indicate that there are aligned grains deep in 
molecular clouds without high mass stars (Ward-Thompson et al. 2000). 

In what follows we discuss grain alignment by radiative torque in molecular clouds. In §2, we calculate efficiency 
of radiative torque in a molecular cloud and minimum aligned grain sizes as a function of visual extinction in the 
cloud. In §3, we calculate polarized far-infrared/submillimcter emission from a pre-stcllar core and discuss the relation 
between the degree of polarization and intensity. We give discussion in §4 and conclusion in §5. 

2. RADIATIVE TORQUES 

As we mentioned above, a flow of photons illuminating a grain can be presented as a superposition of left- and right- 
handed photons, while an irregular grain has different cross section of interaction with photons of different handedness. 
As the result of differential extinction, i.e. absorption and scattering, the grain feels a regular torque. Note, that the 

3 Such a rotation corresponds to the minimum of grain energy for a fixed angular momentum. 

4 It is worth noting that the amplitude of wobbling does not decrease as the efficiency of relaxation increases. The coupling between 
the rotational and vibrational degrees of freedom established by the relaxation mechanism acts back to induce wobbling well in accordance 
with Fluctuation-Dissipation Theorem (see Landau & Lifshitz 1951). More discussion of this point is given in Lazarian & Yan (2003). 



Grain Alignment 



3 



key word here is regular. Random torques produced photons emitted and absorbed by a grain were discussed in terms 
of grain alignment by Harwit (1970). They are rather inefficient, however (Purccll & Spitzer 1971) and are more 
important in terms of damping of grain rotation (see Draine & Lazarian 1998). 

Although the physics of grain spin up by radiative torques was for the most part properly understood by Dolginov & 
Mytrophanov (1976), only calculations in DW96 provided a quantitative insight into the process. These calculations 
obtained for test grains using Discrete Dipole Approximation code (Draine & Flatau 1994) showed that both anisotropic 
and isotropic radiation flows can efficiently spin-up grains. While Dolginov & Mytrophanov (1976) did understand 
that grains will not be only spun up, but also aligned by anisotropic radiation, they could not get correctly what would 
be such an alignment. Numerical simulations In Draine & Weingartner (1997) reveal complex dynamics of grains 
and revealed that in most cases the grains tend to get aligned with long axes perpendicular to magnetic field, even if 
paramagnetic relaxation is absent. Although the nature of this alignment in the absence of analytical calculations still 
remains unclear and some features of grain internal dynamics (that were discovered later!) are missing in the model 
studied (see an attempt in this direction in Weingartner & Draine 2003), it is very plausible that radiative torques 
can provide the alignment that corresponds to observations. Appealing to polarimetric data available one can claim 
that observations do not give us any indications that anisotropic radiation provides alignment that either has wrong 
sign or depends on the angle between magnetic field and anisotropy direction. This would be the case, however, if 
the dynamics of interstellar grains were different from the assumed one. For the rest of the paper we assume that the 
radiative torques do align grains with long axes perpendicular to magnetic field and will concentrate therefore only on 
the magnitude of radiative torques. 

While calculations in DW96 were limited by the interstellar grains, we study radiative alignment of grains of larger 
sizes. Such grains are known to be present in molecular clouds. In addition, unlike DW96, here we are interested in 
the alignment of grains by attenuated and reddened interstellar light that enters into a cloud from outside. 



2.1. Method 



We use the DDSCAT software package (astro-ph/0309069 DW96) to calculate radiative torque on grain particles. 
We use the same grain shape as in DW96, which is an asymmetric assembly of 13 identical cubes. The grain is subject 
to radiative torque because of its irregular shape. We use the refractive index of astronomical silicate (Draine & Lee 
1984; Draine 1985; Loar & Draine 1993; see also Weingartner & Draine 2001). 

In our calculations, the incoming radiation is parallel to the principal axis ai of the grain. Therefore the target 
orientation angle 0, the angle between the incident radiation and the grains primary axis ai (see DW96), is zero. 
Therefore in our calculations the radiative torque, IVadi is parallel to cl\ and jrVadl — |r ra( ^ • &i|. 

For a given wavelength and a grain size, the DDSCAT package returns the torque efficiency vector Qr: 

Qr^ 2 Tra \ /9 , (1) 

where T rad is the radiative torque, a e ff = (SV/At:) 1 ^ 3 the effective target radius, V the volume of the target, u rad 
the energy density of the incident radiation, and A the wavelength. When we consider a radiation field with the mean 
intensity J\, the radiation torque becomes 

dA(4^J A /c)A Qr! (2) 



l eff 



where we used u ra d = 47r J\/c. When we perform the summation over A-axis in natural logarithmic scale, the summa- 
tion becomes 

T rad = 2.303A(log 10 X)(a 2 eff /2c) £(47rJ Aii )A?Qr,i, (3) 

i 

where we used A(log 10 A) = 2.303 dX/X. Figure Q] shows that the value of Ai|Qr,i| (= Ai|Qr,i • ai|) and Af | Qr,* | 
(= AflQr.i • ai|) as a function of A/a for large grains. The quantity Aj|Qr,i| is useful for integration in equation J5J 
and Af|Qr.i| for integration in equation J3J. 

Mathis, Mezger, & Panagia (1983) showed that the average interstellar radiation field (ISRF) in the solar neigh- 
borhood consists of a small UV component plus three blackbody components with T = 3000, 4000, and 7500 K. The 
blackbody components are given by 

"''^^"T ^MTi ' < 4 » 

where W(T = 3000) = 4 x 1Q- 13 , W(T = 4000) = 1.65 x 10~ 13 , W(T = 7500) = 1 x 10~ 14 , k = 1.38 x 10~ 16 , and 
h = 6.63 x 10~ 27 in cgs units. See Figure 1 of Mathis et al. (1983). Radiation field inside a giant molecular cloud 
located at tq = hkpc is also given in Figure 4 of Mathis et al. (1983). They considered a spherical Giant molecular 
cloud that has an isotropic radiation (i.e. ISRF) incident upon the surface of the cloud. They produced the mean 
radiation intensity J\ as a function of the visual extinction Ay measured from the surface of an opaque cloud. We 
calculate radiative torque inside a giant molecular cloud located at tq — hkpc using the radiation field given in Mathis 
et al. (1983). 
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Once we know Qr,i (from DDSCAT) and J\.i (from Mathis et al. 1983), we can obtain the torque from equation 
(0. The gas drag damps grain angular rotation. The gas drag torque is given by 

drag. gas ' | = (2/3)5n_f J -(1.2)(87rm_H-fcT) 1/2 ag // w, (5) 

where tih is the hydrogen number density ran is the mass of a hydrogen atom, S ~ 2, and lo is the angular frequency 
(seeDW96). By equating the radiative torque \r ra d\ (= | ~^rad'3-i\) inequation (0) and the gas drag torque |rdra<?,c/as'^i| 
above, we can obtain the angular velocity of grain rotation around a^: 



LOrad 



■ rad 



(2/3)i5n/f(1.2)(87rm ff fcT) 1 /2 a 4^ yi + T drag ^ as /T drag 



(6) 



where values and definitions of the gas drag time and thermal emission drag time, Td ra g,gas and Tdrag,em respectively, 
are given in DW96 5 . 

The thermal rotation rate lot is the rate at which the rotational kinetic energy of a grain is equal to kT/2: 

= s — ■ (7) 

When a grain rotates much faster than lot, the randomization of a grain by gaseous collisions is reduced. Therefore, 
if a grain rotates superthermally, we expect the grain rotation axis a.i is aligned with magnetic field. From equations 
and (J7J), we have 



LO r ad 
LOT 



LOrad 
LOT 



|r,. 



ad 1 



{2/3)5n H (l.2)(8TTm H kT) 1 / 2 a 



eff , 



1 



5a\ 

192(5 2 \ n H kT 



P a eff 
m H 




7 / dAQ r A(47rJ A /c 



(8) 



(9) 



where 7 is the anisotropy factor of the radiation field. We use 7 = 0.1 for diffuse cloud and 7 = 0.7 for the GMCs in 
accordance with DW96. When the ratio is larger than 1, radiative torque is an efficient mechanism for grain alignment. 
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Fig. 1. — Dependence of A|Qr ■ &l| (left panel) and A 2 |Qr ■ ai| (right panel) on A/a, where A is the wave length and a is the grain size. 
These quantities are useful for estimating which part of the electromagnetic spectrum contributes most to the radiative torque. 



2.2. Results 

We show the results for (Lo ra d/LOT) 2 in Figure|21 The solid line is for the interstellar radiation field (ISRF) in the solar 
vicinity (see Mathis et al. (1983) for details about the radiation field). DW96 used this radiation field and obtained 
the Lo r ad/LOT ratio for three grain sizes (o e // = 0.02,0.05, and 0.2/xm.). In their calculation, they included both 
isotropic and anisotropic components of radiation field. Our calculations arc slightly different. Indeed, we consider 
only anisotropic radiation component the effect of which on alignment is substantially stronger for the grain that 
we use than that of the isotropic component (see Table 4 in DW96). 6 The calculations of radiation anisotropy in a 
turbulent molecular cloud made for us by Tom Bethel show that we do not overestimate 7s. On the contrary, these 
calculations testify that in this paper, if anything, we underestimate the actual values of radiative torques. 

5 Additional processes, e.g. plasma drag, were discussed in Draine & Lazarian (1998). These processes are essential for small grains, but 
less important for large grains that we primary deal with here. 

6 We believe that this is generally true for an ensemble of grains of arbitrary shapes, but more studies are necessary to prove this point. 
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Another simplification is that we consider only anisotropy of radiation only along magnetic field. This is justifiable 
for obtaining a crude estimate, which is the actual goal of our paper. In addition, unlike DW96, we use the UV 
smoothed refractive index of silicon (see Weingartner & Draine 2001). Nevertheless, our result (solid line) agrees with 
that of DW96 within a factor of ~2. 

In Figure we show aligned grain size as a function of visual extinction Ay. We used Figure [3 and assumed that 
grains with u) ra d/u>T > 5 are aligned. For a cloud with n = 10 , 0.2/im grains are aligned even at A v ~ 4. However, for 
a cloud with n = 10 5 , 0.2/im grains are hardly aligned. Grains of ~ l/xm are aligned even at Ay ~ 10 if the density 
does not exceed several times 10 5 cto~ 3 . 

In their classical paper, Mathis, Rumpl, & Nordsieck (1977; hereafter MRN), constructed a model for size distribution 
of dust grains in the diffuse interstellar medium (ISM). This MRN distribution has a sharp upper cutoff at a rnax = 
0.25/im. The MRN model provides a good fit to interstellar extinction and therefore widely used for modeling the 
diffuse ISM. It is expected that at larger optical depths the upper cutoff occurs at larger values. For example, Kim, 
Martin, & Hendry (1994) used the maximum entropy method and obtained a smooth decrease of size distribution 
starting at 0.2/xm. Weingartner & Draine (2001) also obtained an extended distribution beyond the MRN upper 
cutoff. Physically, coagulation of grains happens in denser parts of the interstellar gas (see discussion in Yan & 
Lazarian 2003, Yan, Lazarian & Draine 2004). Therefore it is resonable to assume that grains larger than the the 
usual MRN cutoff are present . Since coagulation is more frequent in dense clouds than the diffuse ISM, we expect to 
see a substantial amount of grains larger than 0.25/xm in dense clouds (see, for example, Clayton & Mathis 1988; Vrba, 
Coyne, & Tapia 1993). If grains of ~ 1/im are abundant in dark clouds, they can emit polarized infrared radiation 
even deep inside the cloud. 
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Fig. 2. — Efficiency of radiative torque. When u) ra( ^/u)T > 1, radiative torque can rotate grains superthermally, which results in grain 
alignment. Different curves represent radiative torque by different radiation fields. The visual extinction Ay is for a giant molecular 
cloud located at 5kpc from the Galactic center. We assume njj = 10 4 cm -3 and T = 20K for the GMC (see Table 6 in DW96 for other 
parameters). For diffuse fSM, we use riu = 30cm -1 and T = WOK (see Table 5 in DW96). 



Fig. 3. — Minimum aligned grain size vs. visual extinction Ay. 
consider 3 different densities. 



We use T = 20K parameters given in DW96. However, note that we 



2.3. Polarization: Rayleigh reduction factor 

In Figure 01 we plot the Rayleigh polarization reduction factor R (p. 328 in Greenberg, 1963; see also Lee & Draine 
1985), which is a measure of imperfect alignment of the grain axes with respect to magnetic field. The conventional 
definition of the factor is R = 1.5((cos 2 j3) — 1/3), where (3 is the angle between the grain angular momentum vector 
and magnetic field. The degree of polarization is reduced when some grains are not perfectly aligned in respect to 
magnetic field. In our case this happens for an ensemble of grains some of which, namely, small ones, are not aligned, 
while the other, namely, the large ones, are perfectly alinged. As the polarization for the range of far infrared A and 
grain sizes a does not depend on those parameters, we can calculate the reduction factor for the entire distribution of 
grains as follows: 

f max C ran n{a)da 

p _ J aligned ™" ( \ (\\ 

Ja min ^ra n n(a)da 

where C ran is the cross-section, n(a) the grain number density, a the grain size, a m i„ the minimum size of grains, 
a m ax the maximum size, and a a Ugned the minimum aligned size, which is given in Figure |3 We assume MRN-type 
power-law grain size distributions: 

n(a)(xa- 3 - 5 (11) 

for a — 0.005/im to a = a max fim. We consider two values for a max : the original MRN cutoff at a max = 0.25/im and 
a larger cutoff at a max = 1/im for the calculation of R. We show the results in Figure Ufa) and^Jb), respectively. 
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For the original MRN distribution (Figure^Ja)), R is smaller than ~ 0.4 for nn > 10 4 cm~ 3 . However, for the larger 
upper cutoff (Figure^b)), R is about ~ 0.2 inside clouds at Ay = 10 when nn ~ 10 5 cm~ 3 . 

We present the results for an opaque giant molecular cloud located at 5kpc from the Galactic center. As we explained 
earlier, we used radiation field given in Mathis et al. (1983). They calculated the radiation field assuming that the 
visual extinction Ay at the center measured from the surface is 200. However, as long as the central visual extinction 
is larger than ~ 15, the radiation field may not be sensitive to the choice of the central Ay (see Flannery, Roberge, & 
Rybicki, 1980). Therefore, our qualitative results obtained here are applicable to various astrophysical objects - from 
dense prestellar cores to giant molecular clouds. 




Fig. 4. — Rayleigh polarization reduction factor (see equation 1101 for the definition in our case), (a) The original MRN distribution 
with cimax = 0.25/im. (b) An extended MRN distribution with a max = 1.0/xm. 



3. POLARIZED EMISSION FROM A DARK CORE 

In the previous section, we showed that larger grains can rotate super-thermally even at Ay = 10 in giant molecular 
clouds. In this section, we apply the result to dense pre-stellar cores. As we noted at the end of the last section, we 
obtained the results in the previous section using the radiation field suitable for giant molecular clouds. Therefore, it 
is questionable whether or not we can directly apply the results in the previous section to prestellar cores. However, 
judging from Flannery et al. (1980) calculation, we expect that the direct application is ill-justified only near the very 
center of the cores. 

3.1. Method 

In this section we calculate polarized emission from a dark pre-stellar core. We assume a simple spherically symmetric 
density distribution and a constant temperature (T ~ 20K). We take density profile of logatropic sphere (Lizano & 
Shu 1989; McLaughlin & Pudritz 1996), which has a finite central density and a p cx r -1 envelope. The logatropic 
sphere is supported by turbulent pressure (and isothermal gas pressure in its original form). The turbulent pressure 
represents nonthermal velocity dispersion observed in clouds. The density profile of a logatropic sphere is different 
from the critically stable isothermal Bonner-Ebert sphere (Bonner 1956; Ebert 1955), which has a p cx r -2 envelope. 
Although some observations (e.g. Alves et al. 1998; Lada, Alves, & Lada 1999; Johnstone & Bally 1999) support the 
p cx r~ 2 profile, other observations (e.g. van der Tak et al. 2000; Colome, di Francesco, & Harvey 1996; Henning et 
al. 1998) support the other profile. For simplicity, we use 

, . f constant if r < ro/4.7 
p(r) cx < _i . (12) 

ry ' I r otherwise, v ' 

where ro is a parameter in our calculations (see McLaughlin & Pudritz 1996 for its physical meaning) and we set 
the central number density nn.c to 3 x 10 5 cm -3 . This distribution truncates at r ~ 24rrj. We take a magnetic field 
from our earlier direct 3-dimensional numerical simulation (see Cho & Lazarian 2003). Numerical resolution is 216 3 
and the average Mach number is ~ 7. The magnetic field has both uniform and random components. The strength 
of the mean field is about 2 times stronger than the fluctuating magnetic field. We assume that the uniform field is 
perpendicular to the line of sight of the observer. 

We assume that the visual extinction A v at the center measured from the surface is ~ 10. This means that the total 
column density through the center is about Njj ~ 3.7 x 10 22 cm~ 2 . The size of the cloud corresponds to ~ 0.02pc. 
This cloud is similar to, for example, L183 (see Crutcher et al. 2004). 

We assume an MRN type grain size distribution, n(a) cx a~ 3 ' 5 , from a = 0.005/im to a = a max . Unlike the original 
MRN distribution, where a max = 0.25/im, we use a max of up to 2/im. We assume that grains are oblate spheroid with 
the axial ratio of ~ 1.2, which is smaller than the value used by Padoan et al. (2001). 
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We follow the method somewhat modified from the one in Fiege & Pudritz (2000) to compute the polarization maps. 
Here we briefly describe the procedures. It is natural to assume (see Fiege & Pudritz 2000) that one can ignore the 
effects of absorption and scattering when dealing with submillimetcr wavelengths. Therefore the polarization submm 
range is due to pure emission. The Stokes parameters are given by 



I cc C r 



, Cp i R 
pels - — / p 



Q oc C P oiRq, 
U oc CpoiRu, 

cos 2 7 1 

~ 3 



ds 



where 



Cpoi — C± C|| , 

C ran = (2Cj_ + C||)/3, 

poos 2 tp cos 2 jds, 



psm ip cos 2 7<is, 



(13) 
(14) 

(15) 



(16) 
(17) 

(18) 
(19) 



R is the polarization reduction factor, tp the angle between the projection of the local B on the plane of the sky and 
north, and 7 is the angle between the local B and the plane of the sky. As explained earlier, the factor R is the n 
reduction factor due to imperfect grain alignment. Note that a a Ugned is a function of both nu and Ay. Figure shows 
how aaUgned is related to nn and Ay. Based on the figure we assume that 



a a iigned = [log 10 (n H )} 3 (A v + 5)/2800 /im. 



(20) 



The error of this fitting formula is around ~ 10%. Note that this fitting formula does not have any physical background. 
From Q, U, and /, we can obtain polarization angle x an( i the degree of polarization as follows: 

tan2x = u/g, (21) 
P= V^. (22) 




Fig. 5. — Polarization map and p-I scatter diagram for the original MRN distribution (i.e. a max = 0.25). We use a logatropic sphere 
for density, which has a p oc r _1 envelop, (a) From the center to the boundary, contours represent 90%, 50%, and 10% of the maximum 
intensity, (b) The scatter diagram follows p oc 7 1 near the center. 
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Fig. 6. — The slope gets flatter as the the upper cutoff, a max , increases. This is because larger cutoff means more aligned grains near 
the cloud center and, therefore, higher polarization intensity, (a) a max = 0.25/xm. (b) a ma x = 0.5/^m. (c) a max = 2.0/xm. 
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Fig. 7. — Polarization map and p-I scatter diagram for the original MRN distribution (i.e. a ma x = 0.25). We use the isothermal 
Bonncr-Ebert sphere for density, which has a p oc r~ 2 envelop, (a) Polarization map. Contours represent 90%, 50%, and 10% of the 
maximum intensity, (b) p-I relation. 



3.2. Simulated map and p-I relation 

In Fig. |Sfa), we plot the polarization map for the original MRN distribution with a max = 0.25/im. We clearly 
observe the depolarization effect, namely, the degree of polarization decreases toward the cloud center. On the other 
hand, the intensity of dust IR emission is strongest at the center. Three contours represent 90%, 50%, and 10% of the 
maximum intensity. In Figure [Sfb), we plot the degree of polarization (p) vs. intensity (J). 

Many observations show that the degree of polarization (p) decreases as intensity (I) increases. The relation is 
usually fitted by a power law. However, the exact power law index varies. Matthews & Wilson (2000) reported 
p oc i~ 7 for the OMC-3 region of the Orion A. On the other hand, Matthews & Wilson (2002) obtained p oc J -0 ' 8 
for the dense cores in Barnard 1. For other dense cores, Henning et al. (2001) reported p oc J -0 - 6 , Lai et al. (2002) 
p oc /~ - 8 , and Crutcher et al. (2004) p oc J -12 . See Figure 1 of Goncalves, Galli, & Walmsley (2004) for illustration. 

We claim that the power law index is sensitive to the value of a max - In Figure H3 we show the change of slopes as a 
function of the upper cutoff a max . The flattening of the slope can be understood as follows. When a max gets larger, 
large grains become relatively more abundant. Since larger grains are aligned even near the cloud center, we expect 
higher degree of polarization near the center. Therefore, the depolarization effect becomes less pronounced and the 
slope gets flatter. We leave for further studies to establish whether or not the actual slope of the curve can be used to 
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constrain the grain size distribution. 

For the sake of completeness, we also calculate the polarization map and p-I diagram for the isothermal Bonncr- 
Ebert sphere. We take the same central density and other parameters as in the logatropic sphere. In Figure Efa), we 
observe more pronounced depolarization effect, which is because of steeper density gradient. The p-I scatter diagram 
in Figure Efb) reflects this. 

4. DISCUSSION 

Our calculations show a substantial increase of radiative torque efficiency as the grain size grows. For the power-law 
distribution of grain sizes we have shown that the p-I relation is sensitive to a max and the density profile of clouds. 
Power-law distribution as an approximation for an actual grain size distirbution was used only for illustrative purpose. 
For instance, in Weingartner & Draine (2001) the grain distribution is approximated by a power law up to a max and 
a exponential tail of grains larger than a max . It is clear from our calculations that if grains with the size of a max are 
aligned, the grains within the exponential tail should also be aligned. In fact, the answer to a very important question 
whether far infrared polarization reflects the structure of magnetic field at high optical depths does not depends on 
the details of the assumed distribution of grains. It is only essential that a substantial percentage of dust mass be 
in sufficiently large grains. If we want to predict a polarization spectrum (see Hildebrand et al. 2000) then the 
size distribution of grains would matter. For instance, grains of different sizes can have different temperatures and 
contribute to either polarized or not polarized parts of radiation. For clouds with embedded stars, the polarization 
spectrum would depend on the distribution of stars as well. 

In practical terms our major result is that large grains must be aligned even at high optical depths. This make 
sub-millimeter polarimetry (see Novak et al. 2003) a useful tool for studies of magnetic fields through the entire 
process of star formation. An earlier understanding reflected in, for instance, Lazarian, Goodman & Myers (1997) 
was that embedded stars are essential for revealing structure of magnetic fields at large optical depth. This meant 
that polarimetry might not be able to reveal the role of magnetic fields at the initial stages of star formation. It 
worth mentioning, that observational evidence that grains are aligned in the conditions when the radiation field is 
substantially reduced have been recently claimed to be a major challenge to grain alignment theory (see Goncalves, 
Galli & Walmsley 2005). 

How can we explain that the optical and near infrared polarimetry does not detect an appreciable polarization 
sygnal originating at high optical depth? We believe that this stems from the fact that the optical and near infrared 
extinction is biased towards small grains which are not aligned. Qualitatively the nature of the bias can be understood 
if one recalls that for A < 1ixa c the efficiency of the grains in producing polarized sygnal drops. At the same time the 
grains with a > a c continue to extinct light. Therefore if a substantial number of grains are larger than a c the dichroic 
properties of the medium in terms of the tranmited light are affected only by grains from a a u g ned given by eq. (|20|l to 
~ a c . At the same time for grain emission at A ^> 2ira the degree of polarization is determined by the grains from 

^aligned to &max- 

To illustrate the situation when both near infrared and far infrared polarimetry would show similar results consider 
a case of the grain alignment at the interface of the diffuse-dense cloud described in Whittet et al. (2001). For 
the range of near-infrared measurements from 0.35 /um to 2.2 /Ltm the the optical cross section of grains less than 
0.25 /mi is still proportional to a 3 . For the case of the Taurus Dark Clouds Whittet et al. (2001) showed that for 
low optical extinctions, i.e. < A v < 3 the ratio of the total to selective extinction stays similar to the value of it 
in the diffuse gas, i.e. R v ~ 3, while the wavelength of maximal polarization \ ma x that enters Serkowski law (i.e. 
Px/Pmax = exp [— Kln 2 {X max /X)\ ), increases. Whittet et al. (2001) interpreted this as the result of the of size- 
dependent variations in grain alignment. Lazarian (2003) explained this as the consequence of the radiative torques 
which fail to align small grains at higher optical depths. Our results here support this conclusion. Indeed, if we adopt 
the grain-size distribution with the original cut-off corresponding to the diffuse medium, i.e. a ma x = 0.25 fim, using 
our Figure 3 and ea. (jlUf) we get Raylegh reduction factor (or effective alignment) around ten percent at A v = 3, which 
is a substantial reduction from the value w f for A v of 1. 

At high optical depth the grain-size distributions are rather uncertain. To illustrate the importance of grain size 
growth for alignment, let us use the distribution in WD0I corresponding to R v — 5.5 for A v of 10 and nn = 10 cm -3 . 
According to Fig. 3 only grains with a > 0.6 /j,m are aligned. According to WD01 the favored distribution of silicate 
grains is cut-off at a smaller grain size. Therefore they are not aligned. The carbonaceous grains have a distribution 
with a cut-off at ~ 1 /im. As the result we expect R carb ~ 0.4, which is larger that the value of effective alignment for 
A v of 3 in the previous example 8 . If we use another model of WD01 corresponding to R v = 5.5 that has MRN-type 
distribution of carbonaceous grains up to size a ~ 10 /jm, then the R carb gets close to 0.8! With all these uncertainties 
we clearly see that far-infrared polarimetry can get insight into the magnetic field topology at large optical depths. In 
fact, we believe that far-infrared polarimetry allows consistency checks for the models of grain-size distributions. 

All these results are valid for clouds without embedded massive stars. The radiation field is being enhanced in 
the clouds and therefore we expect more aligned grains. If grain size distribution stays the same as in dark clouds, 
we expect to have high degrees of far-infrared polarization but still relatively little polarization in terms of optical 

7 It is important to realize that the increase of the upper cutoff for the grain size happens partially due to coagulation of smaller grains. 
Therefore this could be achieved without mantle growth. Naturally, the models with larger grains, i.e. WD01 do not violate the dust-to-gas 
ratio. 

8 We expect to have polarization of the level of Ri 3 percent for this case in emission. 
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and near-infrared polarimetry. The details of this picture can be tested using polarization spectrum technique in 
Hildcbrand et al (2000). 

We would claim that establishing why some grains are not aligned are as important as determining why other grains 
are aligned. These are two inseparable parts of the grain alignment problem that must be solved to make aligned 
grains a reliable technique for magnetic field study. On the basis of the present work we believe that we can account 
for the polarization arising from dust in dark clouds. Thermal nipping in the presence of the nuclear spin relaxation 
described in Lazarian & Draine (1999b) accounts why small grains are not aligned by Purcell's torques. Therefore we 
believe that we have a qualitative agreement between the theory and observations (cf. Goncalves et al. 2004). We 
provide a qualitative comparison of the theory and observations in another paper, where we do radiative transfer in a 
model of a fractal molecular cloud. 

Can the alignment be higher than we predict? Yes, we dealt only with radiative torques. In fact, calculations in 
Lazarian & Draine (1999b) show that for sufficiently large grains thermal flipping is not important. As the result 
such grains are not thermally trapped and can rotate fast in accordance with Purcell's original predictions. Naturally, 
these grains will be aligned paramagnctically. The requirement for Purcell's torques to work in dark clouds is for a 
residual concentration (a fraction of a percent) of atomic hydrogen to be present or for the grains to have temperatures 
different from gas. In addition, MHD turbulence can move grains mostly perpendicular to magnetic field lines and 
align them (Yan & Lazarian 2003). All these mechanisms are likely to act in unison increasing the alignment of grains 
with longer axes perpendicular to magnetic field lines. 

Our calculations have been motivated by grain alignment in molecular clouds. Large grains are known to be present 
in accretion disks around stars, e.g. protoplanetary disks. Our work is suggestive that such grains should be aligned 
and therefore reflect the structure of magnetic field in disks. As magnetic fields are believed to play an important role 
in accretion, the importance of this is difficult to overestimate. 

The limitation of our calculation is that we used a magnetic field from a homogeneous MHD simulation without 
self-gravity. In reality, the magnetic field near dense clouds can be very different from the one we used here. Recent 
calculation by Goncalves et al. (2004) shows that hour glass type magnetic field combined with a torus-like density 
profile can cause a depolarized emission from the cloud center. The reduction factor is around 2. Therefore, we expect 
further reduction of polarization when we use a more realistic magnetic field. 

It is also worth mentioning that the radiation fields given in Mathis et al. (1983) is based on the assumption that 
the cloud is spherical and uniform. Real molecular clouds are likely to be inhomogeneous and, possibly, hierarchically 
clumpy. As a result, the radiation has more chances to penetrate deep within molecular clouds (see Mathis, Whitney, 
& Wood, 2002) to allow grains to be aligned at much higher A v . Elsewhere we have obtain the radiation field inside 
inhomogeneous clouds using direct numerical technique similar to the one in Bethell et al. (2004) and intend to improve 
our present work by combining our results here and a more realistic cloud model with realistic radiation field. 

The ability to trace magnetic fields inside molecular clouds is difficult to overestimate. Using Chandrasekhar & 
Fermi (1953) technique one can infer magnetic fields strength both in the cloud and cloud envelope to test whether 
star formation takes place in sub or supercritical regimes (see Crutcher 2004) . The magnetic connection between clouds 
and cloud cores is also essential for understanding the processes star formation. Does magnetic reconnection plays 
important role for removing the magnetic flux from molecular clouds? This can be answered by studies of magnetic 
field topology. In fact, the study of magnetic topology inside molecular clouds could test different models of magnetic 
reconnection, e.g. those discussed in Shay et al. (1998) and Lazarian, Vishniac & Cho (2004). In addition, polarimetry 
studies of magnetic field structure can bring important insight into the structure of MHD turbulence insight molecular 
clouds (sec review by Lazarian & Yan 2005 and references therein). 

Our finding confirm that the present day understanding of grain alignment can account for all the observational data 
currently available. This makes us believe that polarization arising from aligned grains has become a tool based on 
solid theoretical foundations. The latter is important not only for molecular cloud studies but for many other studies, 
e.g. those of comet, circumstellar polarimetry (see Lazarian 2003 and references therein) as well as for predicting 
and entcrpolating to other wavelength the polarized foreground contribution from dust (see review by Lazarian & 
Finkbciner 2004 and references therein). 

5. SUMMARY 

We have studied the efficiency of grain alignment by radiative torque in optically thick clouds. We have estabilishcd 
that the efficiency of radiative torques is a steep function of the grain size. As the result, even deep inside giant molecular 
clouds (Ay < 10), large grains can be aligned by radiative torque. This means that far-infrared/submillimctcr 
polarimetry can reliably reflect the structure of magnetic field deep inside molecular clouds. Our results show that the 
grain size is important for determining the relation between the degree of polarization and intensity from molecular 
cloud dust. 

This work is supported by NSF grant AST-0243156 and the NSF Center for Magnetic Self-Organization in Laboratory 
and Astrophysical Plasmas. This work utilized CITA supercomputing facilities during its early stages. We thank Bruce 
Draine, Dick Crutcher, Roger Hildebrand, John Mathis, and Giles Novak for useful discussions. 
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